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^ ■ Abstract 

O 

pq ■ Exchanging light pulses to perform accurate space-time positioning is a paradigmatic issue of 

\^ . 

Oj. physics. It is ultimately limited by the quantum nature of light, which introduces fluctuations in 

< \ ..e opt.. _.e.en. and .».. to t.e ..o.>ed St..d„. Qua.t... ... (SQL) Q.^l We 

propose a new scheme combining homodyne detection and mode-locked femtosecond lasers that 
f^ . lead to a new SQL in time transfer, potentially reaching the yoctosecond range (10~^^ — 10~^^ s). 

^ , We prove that no other measurement strategy can lead to better sensitivity with shot noise limited 

^ ' light. We then demonstrate that this already very low SQL can be overcome using appropriately 



o^: 



o 

O 



% 



multimode squeezed light. Benefitting from the large number of photons used in the experiment 
and from the optimal choice of both the detection strategy and of the quantum resource, the 
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f— N ' proposed scheme represents a significant potential improvement in space-time positioning. 
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Accurate spacetime positioning has become a crucial issue for future space experiments 
which require increasing resolution over large distances (see for example [J]). The position in 
space (by ranging to a reference) or time (by clock synchronization with a reference) between 
two observers A and B may be achieved through the Einstein protocol which consists to 



B 



repeatedly exchange light pulses [5|]. From a fundamental point of view, this procedure is 
at the root of Einstein's concept of space and time. From a more practical point of view, 
it permits to distribute the time standard over the whole earth and to precisely know the 
relative position of different satellites in space. 

The basic principle relies on the property that, in the absence of dispersion, each pulse 
carries along its propagation a mean light cone variable u = t± x/c which remains constant 
so that the measurement of the time of arrival of each pulse allows either a determination 
of distance or clock synchronization. The generic situation considered in this paper is the 
following (see figure ([1])) : observer A regularly emits light pulses at a rate synchronized to 
its local clock; B receives these pulses and determine their times of arrival by measuring the 
difference between the arrival times of the incoming light pulses and light pulses delivered 
by a source located in B and synchronized to a reference clock in B. The accuracy of this 
measurement relies therefore on the precision of the clocks in A and B and on the sensitivity 
of the determination of the delay between two light pulses, that we will show how to optimize 
in the present paper. 

Such a delay can be measured by at least two ways: the first one consists in measuring 
the arrival time of the maximum of the pulse envelope. We will refer to this procedure as 
a time-of-fiight (tof) measurement. The second method consists in using the information 
contained in the phase of the electric field oscillation by making an interference pattern 
between the pulses arriving from A and a Local Oscillator (LO) derived from the local clock 
in B. This pattern will give the desired information if the phase of the pulse coming from A 
and the phase of the LO in B are locked to their respective local clocks. This method will 
be referred to as a phase (ph) measurement. 

These measurement schemes suffer from quantum limits associated with the quantum 
nature of light [ll]. For a coherent light pulse of central frequency coq and frequency spread 
Auj, quantum fluctuations lead to the so called Standard Quantum Limit (SQL) of ranging 



for either time-of-flight [2] or phase [y, |7| measurements. Those expressions are given by 

{Auf°'= ^ , {Au)ll^ = ^. (1) 

Where N is the total number of photons measured in the experiment during the detection 
time. Let us briefly discuss those two SQL. First, it is clear on these expressions that the 
SQL can be as small as needed if one can use intense enough light, but there are obvious 
practical limitations to the energy carried by the light pulses. In contrast, isolated photons 
give rise to very low photon fluxes, and the corresponding SQL is very quickly a limitation 
of experimental protocols using photon counting techniques. The expressions also show that 
optical frequencies lead to much smaller SQL than microwave frequencies because of a larger 
ujq and Auj. Finally, as ujq > Au, the phase method has a better ultimate sensitivity than 
the time-of-flight technique but requires highly spatially and temporally coherent sources. 

For the time being, the resolution in time transfer is limited by classical technical noises 
so that the previous SQL are not yet a limitation in time transfer. Nevertheless, with the 
recent developments in stabilization of frequency combs referenced to optical standard, it 
is getting closer and closer to these quantum limits [8i|. Both for a fundamental point of 
view and for future experiments, it is therefore necessary to compute the ultimate sensitivity 
in time transfer with mode locked femtosecond laser since the latter combine both a time- 
of-fiight information in their enveloppe, and a well stabilized phase information inside the 
enveloppe. 

In order to compute the SQL in timing involving mode locked femtosecond lasers, we 
begin by writing the positive frequency electric field operator ELJ emitted by A in the 
absence of any perturbations, as a decomposition in temporal modes : 



A^oV(") = ^ 5Z ^nVn{u) , S = iJ^^ , (2) 

n 

where T is the measurement time. The orthonormal temporal modes f„(u) will be written 
as a (complex) time- varying amplitude gn{u) multiplied by a propagation phase factor of 
the form e~*'^°" : 

Vn{u)=gn{u)e-''^°^. (3) 

The annihilation operator corresponding to those modes are noted d„. Without any loss 
of generality, we can appropriately choose the mode basis such that the mean value of the 



electric field operator ELMu) is proportional to vq, namely (ELMu)) = EyNe^^ vq{u), with 
A^ the mean number of photon and 9 a global phase. 

Now, any variation An of the mean light cone variable, caused for example by a distance 
change between A and B, leads to a modification of the field received in B which reads 
&^\u) = ELMu — Au) (see figure ([I])). The temporal mode corresponding to this field can 
be decomposed as follows if the perturbation Au is small : 



Vo{u — Au) ^ Vo[u) — Au 



= Wo(m)H wi{u) . (4) 

u=0 «o 



du 

The constant Uq ensures the normalization of the new mode Wi{u). The latter one will be 
called the timing mode because it carries the timing signal Am. For pulses of frequency 



spread Aa; [39|], Uq is given by Uq = 1/^/uJq + Auj'^ and the expression of the timing mode is 



^1*^") = / o , 1 [io^Mu) + Mu)) , a = -^ . (5) 

a is roughly equal to the number of field oscillations within the pulse, which can be as 
small as a few units for femtosecond pulses. The timing mode Wi{u) contains two terms: 
the first one, namely wo{u), gives a contribution to the timing signal via a phase change 
(interferometric method of ranging). The second one, namely vi{u), is normalized and 
orthogonal to vq so that it will be taken as the second mode of the basis {vn)n- It reads : 

.,H = --^%Me--V (6) 

^ ^ Aiu du ^ ^ 

This mode gives a contribution to the timing signal via a time shift of the pulse enveloppe 
(time-of-fiight technique). The latter mode is represented in the figure ([2]) and is the tem- 
poral analog of the spatial TEMqi gaussian mode when the emitted pulses are gaussian. 

The timing signal Au can be retrieved by projecting Vq{u — Au) on the timing mode 
Wi{u). This can be done using the balanced homodyne detection scheme represented in 
figure (I2l) where the input pulses are mixed with a Local Oscillator (LO) put in the timing 
mode wi 9], so that {EIq{u)) = Sy/Ni^oe^^^°wi{u), with A^lo the mean number of photon 
in the LO field and ^lo its phase. Denoting (&„)„ the annihilation operators for the LO, the 
homodyne signal D reads : 

D=\S\'Y,{^iK + bia^). (7) 



The mean signal of the balanced homodyne detection when a timing offset Am is present, 
then reads : 



{D) = 2\SY^NN] 



LO 



A7/, 01 

cos(0 - 0Lo) + r-o T- sin(^ - ^lo) 



(8) 



We assume from now on that, as usual, the LO is much more intense than the input 
field. The general case can be treated without difficulty. In this situation, the variance of 
the balanced homodyne signal, taken for Au = 0, is given by : 



\£rN, 
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where a\ and a^ are the variances of the quadrature operators Pq (phase operator of mode 
Vq) and Qi (amphtude operator of mode Vi) of the input field : 

Po = i (ale'^^° - aoe"*^^o) and Qi = ale'^'^^ + aie''^"^" . (10) 

The Standard Quantum Limit (SQL) is then obtained as the smallest Am that can be 
measured using shot noise limited coherent light (cr|, =0"^ =1); assuming a signal to noise 

^0 Wl 

ratio equal to one {{D) = ajj). It is obtained for 6 = ^lo and is given by : 

The expression ( TTTj) is one the main results of this paper and gives a new SQL in timing. The 
latter is lower than both the SQL in time-of-flight and phase measurements (see equation 
([1])), which obviously are special cases of our scheme when the LO is either in the ivo 
or vi mode. This means that the proposed balanced homodyne detection scheme has a 
better sensitivity than existing schemes based on either time-of-flight or interferometric 
measurement. The improvement comes from the fact that coherent pulses, in addition 
to their phase, carries a time of flight information in their time varying enveloppe. Both 
informations are read by the balanced homodyne detection if the LO is shaped in the mode 
Wi{u). Let us stress that such optimized measurements have already been successfully 
employed for pure phase measurement [10| and in the spatial domain to measure transverse 



beam displacement and tilt 11|. 



For a P = 10 mW laser with A ~ 810 nm and a 10 fs pulse duration, the SQL is equal to 
(Am)sql = 2 X 10^^'^ s, i.e. a noise level of 2 x 10^^^ s/v^Hz (20 yoctoseconds for one second 
integration time). 



A natural question is to know whether it is possible to reach still better sensitivity on the 
same beam but by using another measurement strategy. An answer can be provided in the 
context of information theory with the help of the Cramer- Rao bound [12] , which gives the 
smallest measurable delay Am that can be achieved in the presence of a given distribution 
of noise. This bound has the property of being independent of the measurement strategy 
and depends only on the noise of the incoming signal. A calculation of the Cramer-Rao 
bound, analogous to the one detailed in [l3|, [ij] proves that using coherent light, this bound 
is precisely equal to the expression (TTT!) of (A'u)sql- We are therefore sure that no other 
measurement scheme will reach a better accuracy than the introduced balanced homodyne 
detection and in this sense this scheme is said to be efficient. 

Obviously the SQL flTTl) is the fundamental limit when one restricts oneself to the use of 
classical states of light and coherent states, as proven with the previous standard Cramer- 
Rao bound. Nevertheless, it is well known that it can be beaten using quantum resources 3|, 
llSl . [iGj. For example, the improvement of the sensitivity in interferometric measurements 
using squeezed light has been proposed p, LZl] , observed experimentally 171, [l8|, [l9|] , and will 



be certainly practically implemented in the future generations of interferometric detectors 



of gravitational waves 20]. The use of an entangled photon source to improve time-of-flight 



ranging measurements in the photon-counting regime has been also proposed [2, l2l| and 
experimentally demonstrated [22] at a picosecond level of timing sensitivity. We propose 
here to improve the scheme introduced previously by using appropriately squeezed light. 

Inspection of equation ([9]) immediately shows that in the case of a strong LO the signal 
to noise ratio is increased if the noise of the incoming mode wi is below the shot noise. 
This can be obtained if squeezing of the input field modes Vq and Vi is achieved along the 
quadratures Pq and Qi respectively. This therefore requires to first squeeze the phase of 
the input field and mix it with a squeezed vacuum mode Vi [23J, using procedures already 
demonstrated in the spatial domain 91]. If we assume that the squeezing coefficient is equal 
for the two states, namely a-^ = u-^ = e~^ (r > 1 being the squeezing parameter), then 
the new minimum measurable value of Am is given by : 

^^""^^^-^'-^ = vf74+A^ "'" ■ ^^^^ 

This minimum resolvable Am is thus reduced below the SQL (TTT]) by the factor e^. Note 
that the expression for the general case of different squeezing along Xq and Xi, as well as a 



LO not supposed strong, can be obtained straightforwardly from the equations given in the 
paper. 

Using the best present technology, the noise reduction factor can reach 10 dB 
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25|,i.e. 



a factor of 10 improvement, even at low noise frequencies. The advantage of squeezing over 
the other proposed quantum techniques such as entanglement is that it can be used together 
with an intense beam for which the SQL is already very low. In addition the squeezed beam 
travels along with the signal beam and therefore they both share commun noises. This is 
not the case for protocol using non local quantum correlations (such as most of the protocols 
based on entanglement) where the quantum correlations are spread over a large region of 
space and submitted to differential noise effects. The main drawback of squeezing is its 
sensitivity to losses in the optical system and the detectors. This means that the technique 
could be used in situations where light propagates in vacuum, for example between satellites 
in flying formation. 

An experimental implementation of the scheme with the aim at reaching the SQL and then 
observe the quantum improvement suffers different technical challenges. First of all, reaching 
a timing precision in the yoctosecond regime requires very stable laser repetition rate and 
phase stabilization. This can be eventually be achieved with mode-locked femtosecond 
lasers which are already used for absolute and relative ranging in different measurement 
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291]. The dominant source of noise in equation (Q is given by the 



schemes 

noise cr^ of the phase of vq. Self-referencing stabilization using a / — 2/ beat allows to 



keep this noise to a very low level, down to 10 ^rad/vHz at 10^ Hz with state-of-the-art 
stabilization techniques 30, l31|, corresponding to a timing noise of 4 x lO^^^s/vHz at 



10^ Hz. Concerning the repetition rate Trep, the latter can be locked to an optical re: 



and current technology leads to a time jitter noise level of 10 s/v Hz at 10 Hz |32 



erence, 
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Another experimental challenge is to produce the squeezed temporal mode wi. Indeed, if 
the mode Wi can be obtained with presently available commercial mode shapers 35|, the 
squeezing is much more challenging, but can in principle be obtained by propagation through 
a non-linear Kerr medium 36| or more efficiently by using parametric down conversion 
pumped by mode-locked lasers [37|, l38|, or even synchronously pumped OPOs [23 1. 
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FIG. 1: a) General scheme of a one way time transfer. We only consider propagation along the 
X axis with diffraction neglected, b) Spacetime representation in the reference frame of observer 
B (xb = 0). a modification Axa of the position of the observer A leads to a modification An = 
— Axa/c of the light cone variable that is emitted towards B and consequently leads to non regular 
time of arrival in B. 
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FIG. 2: Proposed balanced homodyne scheme to reach optimal detection in ranging measurement. 
The pulses synchronized on the clock in A are measured in B by homodyne detection with pulses 
synchronized on the local clock and in an adequate temporal mode (here is represented only the 
part vi of the LO for clarity). 
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